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Abstract 

P mono-doped and (P, N) co-doped ZnO are investigated by the first-principles 
calculations. It is found that substitutive P defect forms a deep acceptor level at O site 
(Po) and it behaves as a donor at Zn site (Pzn), while interstitial P (Pi) is amphoteric. 
Under equilibrium conditions, these defects contribute little to the p-type conductivity 
of ZnO samples since the formation energy of Pzn is much lower than that of Pi or Pq 
when Ef is below mid-gap (a prerequisite p-type condition). Zinc vacancies (Vzn) and 
Pzn-2Vzn complex are demonstrated to be shallow acceptors with ionization energies 
around 100 meV, but they are easily compensated by Pzn defect. Fortunately, Pzn-4No 
complexes may have lower formation energy than that of Pzn under Zn rich condition 
by proper choices of P and N sources. In addition, the neutral Pzii-3No passive defects 
may form an impurity band right above the valence band maximum of ZnO as in 
earlier reported (Ga,N) or (Zr,N) doped ZnO. This significantly reduces the acceptor 
level of Pzn-4No complexes, and helps improving the p-type conductivity in ZnO. 
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I. INTRODUCTION 

Zinc oxide (ZnO) is a good candidate for applications in short-wavelength and 
transparent optoelectronic devices for its unique physical properties. Yet difficult 
realization of reproducible low resistivity p-type ZnO hinders its applications. There 
have been many attempts, such as mono-doping, [1-10] co-doping or cluster doping 
technology [11-14] to overcome the /'-type doping bottleneck. Among those attempts, 
donor-acceptor co-doping method is promising in terms of enhancing the acceptor 
concentration and lowering the acceptor level. [13, 14] A record of /(-type electrical 
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resistivity of 0.026 Qcm, Hall mobility of 4.4 cm V" s" , and carrier concentration of 
5.5x10 cm" has been reported with zirconium (Zr) and nitrogen (N) co-doping 
ZnO. [14] Reduction of the transition levels through the effective impurity bands 
induced by donor-acceptor complexes have successfully explained the experiment 
observation in Ga-N[15, 16] and Zr-N[16] co-doped p-type ZnO. Recently, co-doping 
with group V element and N, such as As-N [17] and P-N, [18] have also realized 
p-type ZnO. A detailed explanation of this dual-acceptor doping mechanism is not 
available so far, although it could be viewed as an attempt following the earlier 
research on group III-N and group IV-N co-doping. 

In order to clarify the p-type conduction in (P, N) co-doped ZnO samples, P 
mono-doping and (P, N) co-doping in ZnO are studied by first-principles calculations 
in this work. By comparing the defect formation energies, we find that there are three 
kinds of defects, i.e., Vzn, Pzn-2Vzn, and Pzn-4No complexes may be shallow 
acceptors. Furthermore, it is found that acceptors Vzn and Pzii-2Vzn may not be 
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responsible to the high hole carrier concentrations of about 10^^-10^^ cm"^ [3, 4] as 
observed experimentally due to the compensation effects of Pzn under equilibrium 
conditions, no matter what kind of P or N sources are adopted. However, Pzn-4No 
complex may play a critical role in p-type conductivity in (P, N) co-doped ZnO under 
Zn-rich conditions by using P with low A//p and N with highA/z^j dopant sources. 
Special attentions have also paid to the mechanism of the lowered transition level of 
Pzn-4No complex by passive impurity band. 

II. COMPUTIONAL METHODS 

Our calculations are carried out with the Vienna ab initio simulation package 
(VASP), [19, 20] using the generalized-gradient approximation (GGA) formulism of 
PW91 [21, 22] for the exchange correlation. The electron-ion interactions are 
described by the projector augmented wave method (PAW). [23, 24] Zinc 3d states are 
treated as valence electrons. The energy cutoff for the plane wave expansion is set to 
500 eV. All atoms are fiilly relaxed during the calculation until the 
Hellmann-Feynman forces converge to 10 meV/A. The optimal crystal lattice 
constants of ZnO (a=3.28A, c=5.27 A) are adopted through the calculations, which 
are in excellent agreement with the experimental values (a=0.325 nm, c=0.521 nm). A 
72-atom supercell is used to simulate P mono-doping and (P, N) co-doping in ZnO. 
The Brillouin zones are sampled with gamma centered 3x3x2 ^-point mesh. For 
charged defects, a uniform background charge is added to keep the global charge 
neutrality of the supercells. 
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The defect formation energy, Ml j (a,q), for a supercell containing defect a in 

charge state q, can be computed as: [25-27] 

AHf (a, q) - AE(a, ^) + Z n^AjUf + qE^ , ( 1 ) 

where AE (a, q)^E(a,q)-E (host) +'Zn.ju. (hulk) + q £yj^y^(host) . E(a,q) is the 

total energy for the studied supercell containing defect a in charge state q and 
£'(host)is the total energy of the same supercell without the defect. «, is the number 

of atom i involved in the defect, and q is the number of electrons transferred from the 
supercell to the reservoirs in forming the defect cell. E^ is the electron Fermi energy 
referenced to the valence-band maximum (VBM) of the host, and varies up to the 
experimental bang gap value of 3.37 eV. [28] £"vbm is the VBM energy of the host 
supercell. A//, is the chemical potential of atom i referenced to its elemental 
solid/gas with cohesive energy of ju^ (bulk) , and should not be greater than eV 
generally in order to avoid precipitation of the elemental solid/gas. Meanwhile, to 
maintain a stable compound ZnO and to avoid other possible competing phases, it 
must satisfy the following conditions 
Ay«z„ + ^Mo = ^/ (ZnO) 

J ■ ^Mzn +k-AjUo+m-Aju^+n-Aju^<AHf (Zn .0,P^N„ ) (2) 
Where Zn^ O^^P^N^ stands for possible alloys formed by Zn, O, P and/or N naturally 
except for ZnO. A//^(ZnO) and A//^(Zn^.O^P^N„) represent the formation 
energies of ZnO and the corresponding competing phases. The calculated cohesive 
energies of elements and the calculated formation energies of possible competing 
compounds in (P, N) co-doped ZnO are listed in Table 1 . 
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The defect transition energy level s^iq/q ) is the in Eq.(l), at which the 
formation energy AHj. (a, q) of defect a in charge state q is equal to that of another 

charge 9' of the same defect, i.e., 

(q/q) = [AEia, q) - ^E{a,q)-\l{q - q) (3) 
In this paper, we use the hybrid scheme combining A:-points sampling and 
F-point-only approaches to calculate the transition energy level and the defect 
formation energy. [27, 29] In this scheme, the transition energy level for acceptors 
(g<0) with respect to VBM is given by: 

^(0/^) = [4 (0) - ^vBM (host)] + [E{a, q) - {E{a,Q) - qe'^ (0))]/(-^) (4) 
For donors (q>0), the ionization energy referenced to the conduction-band minimum 
(CBM) is given by: 

< (host) - s(0/q) = (host) - 4 (0)] + [E(a, q) - (£(«,0) - g4 (0))]/g (5) 
where (0) and s]^ (0) are the defect levels at the special A:-points (averaged) and 
at the F-point, respectively; £-yB^^(host) and £-^M(host) are the VBM and CBM 
energy, respectively, of the host at the F-point; and £"^(host) is the calculated 
bandgap at the F-point. Meanwhile, the average electrostatic potential at core area of 
oxygen atom far away from the defect is adopted as a reference to determine the VBM 
alignment. 

The formation energy of a charged defect is then given by: 

AHf(a,q) = AHf(a,0)-q£(0/q) + qE^ (6) 

where AH^^ ia,0) is the formation energy of the charge-neutral defect. 

Although a 72-atom supercell is adopted in the calculations, the defect 
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concentration of the studied system is around 10 -10 cm" , which is significantly 
higher than a typical doping concentration. The fictitious high doping concentration 
may imply unwanted impurity interactions, especially for the charged defects. Here a 
band filling correction is employed to get a more reasonable transition level following 
the approach described in Ref 26. To be consistent with earlier theoretical work on 
p-type doping in ZnO, the image charge correction is not employed in this work. In 
fact, we have checked the transition level of Aszn-SVzn to validate our approach on the 
transition level calculation. Our calculated transition level ^(0/-) is 0.14eV for 
Aszn-2Vzn, which is in good agreement with earlier reported value of 0.1 5eV by 
Limpijumnong et al. [30] 

III. Monodoping of P in ZnO 

Figure 1 shows formation energies of defects related to phosphorous 
mono-doping in ZnO under the two extreme conditions: (a) zinc-rich limit, i.e., 
A//2„ = and AjUq = AH ^{ZaO) ; and (b) oxygen-rich limit, i.e., 
AjUq = OandA//2„ = AH {ZnO) . Here we suppose P2O5 is adopted as the P source. 
The slope corresponds to the charge state q as used in Eq.(l)~Eq.(6). A change in the 
slope indicates transition of the charge state. The transition levels are independent of 
the choice of atomic chemical potentials and thus they are the same in Fig. 1(a) and 
Fig. 1(b). Similar to the reported P[31] and As[30] dopants, interstitial P (Pi) is 
amphoteric as diluted point defects. A substitutional P at an O lattice site (Po) forms a 
deep acceptor level of £"(0/-) = 0.96 eV , which is in excellent agreement with earlier 
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calculation of 0.93 eV. [32] A substitutional P at a Zn lattice site (Pzn) behaves as a 
donor. The formation energy of Pzn is much lower than that of Pi and Pq in most 
situations except that the Fermi level is greater than 2.98 eV above the VBM under 
0-rich condition. This means that none of these three point defects could contribute to 
the p-type conductivity remarkably in ZnO. 

Recent calculations have argued that in ZnO doped with P2O5, the dominant 
acceptors are Vzn defects, while the Pzn-2Vzn complex is energetically more favorable 
for the Zn3P2 source under 0-rich conditions. [31] Yet experiment results indicated 
that P atoms replace the O atoms to form Pq rather than forming Pzn-2Vzn complex as 
acceptors. [1] Here, the formation energy of Pzn-2Vzn complex and Vzn are also 
plotted in Fig.l for comparison. Figure 2 indicates the most stable structure of 
Pzn-2Vzn complex, which is consistent with that reported by Lee et.al. [31] Contrary 
to the reported deep transition level (0.55 eV above the VBM), [31] we find that 
Pzn-2Vzn complex may be a shallow acceptor defect, with the ionization energy of 
£•(0/-) = 0. 12 eV . On the other hand, the formation energy of Pzn-2Vzn complex is 
4.62 eV (Zn-rich) and 2.34 eV (0-rich) with P2O5 source for P source, which is much 
greater than that of Pzn when the Fermi level is close to VBM. This indicates that 
Pzn-2Vzn complex may be fully compensated by Pzn when P2O5 source is adopted. 
The calculated transition levels of Vzn are £"(0/-) = 0.091 eV and 
£(-/2-) = 0.42 eV, indicating that Vzn may act as a shallow acceptor. However, its 
formation energy is relatively high compared with that of Pzn under both Zn-rich and 
0-rich conditions when Ep is near VBM. It also suffers the compensating effect of Pzn, 
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thus unlikely to play important role in />-type conductivity, which requires the Fermi 
level of ZnO below the mid gap. 
IV. (P, N) co-doping in ZnO 

Recently, it was proposed that the introduction of mutually passivated impurity 
bands may successfully overcome the doping asymmetry in ZnO by (Ga, N)[15, 16] 
or (Zr, N)[16] co-doping. Consequently, passive defect complexes might be formed in 
group V, N (such as P, N) co-doped ZnO. ActauUy, there are two noticeable 
configurations of (P-3N) complex (as shown in Fig. 3) when P substituting Zn and 3 
N replacing three nearest neighbor O atoms. One is the three nearest in-plane O atoms 
replaced by N atoms [Fig. 3 (a)], and the other is N atoms substituting the nearest 
out-of-plane O atom together with two of the nearest in-plane O atoms [Fig. 3(b)]. The 
configuration shown in Fig.3(a) is considered as the passive stoichiometric (P-3N) 
complex, since it is energetically favored by 0.039 eV than the latter one. The 
calculated total DOS for the pure ZnO host and the passive stoichiometric (P-3N) 
complex doped ZnO are shown in Fig. 4. Unlike the additional fully occupied 
impurity band in (Ga, N)[15, 16] and (Zr, N)[16] co-doped ZnO, the valence band 
appears to be modified by the passive stoichiometric (P-3N) complex, with its VBM 
shifted up by 0.156eV with respect to that of the pure ZnO host. We have carefully 
checked the partial DOS of the passive stoichiometric complexes in (Ga, N) and (P, N) 
co-doped ZnO to investigate the property of the modified valence band in (P, N) 
co-doped ZnO, and found that the electronic states near the VBM of the doped 
systems is mainly constituted of the p states of N dopants and d states of Zn. The p 
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states of N dopants of the passive stoichiometric complexes in (Ga, N) and (P, N) 
co-doped ZnO, as well as the d states of Zn of the two doped systems and the ZnO 
host are compared in Fig. 5. It can be clearly seen that the upper part of the valence 
band in (P, N) co-doped ZnO has the same origin with (Ga, N) co-doped ZnO, i.e., the 
hybridization of the p states of N dopants and d states of Zn of ZnO host. In other 
words, the passive stoichiometric (P-3N) complex also forms an additional fully 
occupied impurity band, although it is not obviously shown in DOS plot as that in (Ga, 
N)[15, 16] and (Zr, N)[16] co-doped ZnO. With the impurity band above VBM of 
pure ZnO, the acceptor level is expected to decrease when additional N atoms are 
introduced into the (P-3N) co-doped ZnO since electrons may transit from the 
impurity band. 

When an additional N is introduced to the above mentioned two configurations 
of (P-3N) complex, we have investigated nine unequal configurations of P-4N 
complex for the additional N atom replacing the nearest and the second nearest 
neighbor O site of Pzn- We find that the configuration with four N atoms substituting 
the first neighbor O sites of the Pzn is energetically favored, while the other 
configurations are at least 0.64 eV higher in total energy. This clearly indicates that 
the doped P and N tend to form Pzn-4No complexes in ZnO. 

For the dopant source of P2O5 and N2, the calculated formation energy of the 
Pzn-4No complex is shown in Fig. 6. The ionization energy of Pzn-4No is 0.271 eV 
above the VBM of ZnO host according to Eq.(4), yet it is reduced to 0.115 eV from 
the impurity band of the passive (Pzn-3No) complex. Fig.6 also shows that the 
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formation of Pzn-4No in Zn-rich condition is -1.28 eV, which is 18.87 eV lower than 
that of O-rich condition, indicating that p-type conductivity is easier to achieve in 
Zn-rich condition when P2O5 and N2 are used as dopants in ZnO simultaneously. 

The formation energies of three kinds of accepters with low transitional level, 
Pzn-4No, Pzn-2Vzn and Vzn, under two extreme conditions, (a) for Zn-rich and (b) for 
O-rich, are plotted in Fig. 7. The most possible killer of acceptors, Pzn, is also shown 
in Fig.7. We can see that for the three acceptors, Pzii-4No is energetically favorable 
under Zn-rich condition; while Vzn is favored under O-rich condition when P2O5 and 
N2 are served as dopant sources. Yet the concentration of hole carriers is expected not 
to be much for the compensation effect of Pzn donor in both situations for P2O5 and N2 
sources (Solid line in Fig.7). 

We notice that it is rare to conduct the (P, N) codoping under P-rich and N-rich 
condition, i.e, the chemical potentials of P and N are in equilibrium with those of P2O5 
and N2, respectively. Additionally, most experimentally reported p-type ZnO samples 
are obtained with NO source instead of N2. Therefore, there are rooms to further 
investigate the jc-type doping possibility of (P, N) by adjusting the chemical potential 
of the dopants. [33] Since the change of A//p will simultaneously shift the formation 
energies of Pzn-2Vzn and Pzn, the adjustment of the chemical potentials of the dopants 
(i.e., P and N) has no influence to the relative formation energy of Pzn-2Vzn to Pzn 
[Fig. 7(b), dashed line]. This means that although Pzn-2Vzn complex may 
energetically more favorable than Vzn, it may not be accounted for the p-type 
conductivity in ZnO merely through adjusting the chemical potential of P and/or N 
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under 0-rich condition. In other words, adjustment of A//p and Aju^ has no impact 
of pinned Fermi level [Fig. 7(b)] under O-rich condition. 

Considering the formation energy difference between Vzn and Pzn ( ^ = ) is 
about 6.30 eV [solid line in Fig.7(b)], it is very difficult to overcome the energy 
difference by adjusting A//p experimentally. However, the situation may be changed 
under Zn-rich condition. When the chemical potential of P source is reduced and that 
of N source is raised (e.g., by using NO, NO2 sources), the formation energy of 
Pzn-4No will decrease while that of Pzn will increase. This will enhance the p-type 
conductivity in ZnO. For example, when the chemical potential of P source is reduced 
to A//p = -3.0 eV , and that of N source raised to Aju^ = 1.7 eV , the pinned Fermi 
level will decrease to around the transition level of Pzii-4No [ shown with dashed line 
in Fig.7(a)], then the compensation effect of Pzn becomes insignificant and /7-type 
conductivity in ZnO will be greatly enhanced. The formation energy of Pzn-4No could 
even be tuned lower than that of the compensator by further adjustment of Aju^ and 
Aju^ when necessary. From this point of view, we tentatively suggest preparing (P, N) 
co-doped p-typQ ZnO under Zn-rich condition, rather than under O-rich condition. 
IV. SUMMARY 

We have investigated the formation of isolated defects and defect complexes in 
P-mono doped and (P, N) co-doped ZnO samples through the first-principles 
calculations. Although Vzn and Pzn-2Vzn complex show low ionization energy, they 
pay little contribution to p-type conductivity of samples for the strong compensation 
effects of Pzn donors. The Pzn-3No passive defects may form an impurity band right 
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above the VBM due to the hybridization of the p states of N and the d states of ZnO, 
as in (Ga,N) and (Zr,N) doped ZnO systems. Thus the ionization energy of Pzn-4No 
complex is reduced from 0.271 eV to 0.115 eV when electrons are transited from the 
top of impurity band. We also find that Pzii-4No complex may overcome the 
compensation effect of Pzn under Zn-rich condition with a proper choice of P, N 
dopant sources, and attribute to the reported /'-type conductivity in (P,N) co-doped 
ZnO system. 
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Figure captions 

Figure 1 (Color online) The formation energies of possible P mono doped defects, Pzn, 
Pi and Po, together with Pzn-2Vzn complex and Vzn as a function of the Fermi energy 
under Zn-rich (a) and 0-rich (b) condition for P2O5 source. 

Figure 2 (color online) Relaxed structure of the most stable Pzn-2Vzn complex 

Figure 3 (Color online) Structures of (P-3N) complex when P substituting Zn and 3 N 
replacing three nearest neighbor O atoms. Configuration (a) denotes the three nearest 
in-plane O atoms replaced by N atoms, and (b) N atoms substituting the nearest 
out-of-plane O atom together with two of the nearest in-plane O atoms. 

Figure 4 (Color online) Calculated total DOS of the pure ZnO host and (P-3N) 
co-doped ZnO 

Figure 5 (Color online) Comparison of the p states of N dopants of the passive 
stoichiometric complexes in (Ga, N) and (P, N) co-doped ZnO (a) and d states of Zn 
of the two doped systems and the ZnO host (b). 

Figure 6 The formation energy of Pzn-4No complex when P2O5 and N2 are used as 
dopants. The vertical dashed line stands for the impurity band due to the passive 
(Pzn-3No) complex. 
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Figure 7 (Color online) The formation energies of Pzn-4No (blue), Pzn-2Vzn (red) 
and Vzn (magenta) under two extreme conditions, (a) for Zn-rich and (b) for O-rich, 
together with the most possible killer of acceptor, Pzn (Green). The solid lines indicate 
the formation energies for P2O5 and N2 sources, while the dashed lines indicate the 
formation energies with A//p=-3.0eV and Aju^=l.7eV. The vertical short 

dashed line stands for the additional impurity band due to the passive (Pzii-3No) 
complex. 
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Tables 

Table 1 The cohesive energies of elements and the fonnation energy of possible 



competing compounds in (P, N) co-doped ZnO, as well as their space group 



Element 


Space Group 


Cohesive energy (eV) 


Zn 


P63/mmz 


1.123 


0(02) 


PI 


2.889 


N(N2) 


PI 


5.328 


P 


CMCA 


3.542 


Compound 




Formation energy (eV) 


ZnU 


roiJVlC 


-3.525 


ZnOi 


P21/A-3 


-3.435 


NO2 


1 

PI 


-0.959 


NO 


PI 


0.724 


N2O 


PI 


A '^11 

4.211 


N2O3 


P212121 


-1.974 


N2O4 


I2/M-3 


-3.681 


N2O5 


P63/mmc 


-3.752 


Zn3P2 


P42/nmcs 


-1.155 


Zn3N2 


I21/A-3 


0.415 


Zn3(P04)2 


C12/C1 


-30.092 


Zn2P207 


PBCM 


-26.089 


ZnP206 


C12/C2 


-21.676 


P2O5 


FDD2 


-16.580 


P4O6 


P121/M1 


-19.334 


P3N5 


C12/C2 


-2.641 
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Fig. 1 Ren-Yu Tian and Yu-Jun Zhao 
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;. 2 Ren-Yu Tian and Yu-Jun Zhao 
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Fig. 3 Ren-Yu Tian and Yu-Jun Zhao 
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Fig. 5 Ren-Yu Tian and Yu-Jun Zhao 
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Fig. 6 Ren-Yu Tian and Yu-Jun Zhao 
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Fig. 7 Ren-Yu Tian and Yu-Jun Zhao 
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